Geometric phases in superconducting qubits beyond the two-level-approximation 
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Geometric phases, which accompany the evolution of a quantum system and depend only on its 
trajectory in state space, are commonly studied in two-level systems. Here, however, we study the 
adiabatic geometric phase in a weakly anharmonic and strongly driven multi-level system, realised as 
a superconducting transmon-type circuit. We measure the contribution of the second excited state 
to the two-level geometric phase and find good agreement with theory treating higher energy levels 
perturbatively. By changing the evolution time, we confirm the independence of the geometric phase 
of time and explore the validity of the adiabatic approximation at the transition to the non-adiabatic 
regime. 

PACS numbers: 03.65.Vf, 03.67.Lx, 42.50.Pq, 85.25.Cp 
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When a quantum mechanical system evolves under a 
time-dependent Hamiltonian, its wavefunction acquires 
a geometric phase in addition to the dynamic phase. ^ 
If the evolution of the Hamiltonian is cyclic and adi- 
abatic, the geometric phase is termed Berry's phase. ^ 
While the dynamic phase is proportional to the time inte- 
gral of the system energy, Berry's phase simply depends 
on the geometry of the path the Hamiltonian traces out 
in its parameter space and is independent of time and 
energy. It has been observed in a variety of systems, 
ranging from photons in optical fibers^ and nuclear mag- 
netic resonance (NMR)^ to superconducting circuits^, 
superconducting charge pumps'" and electronic harmonic 
oscillators.^ The purported resilience of the geometric 
phase against certain types of noise, studied both in the- 
or exp er imeni! ^ I ^ ^ I ^ has raised interest for im- 

plementing geometric gates in quantum information pro- 
cessing. Recently, the geometric aspects of multi-level 
systems have attracted increased attention. The geomet- 
ric phase has been observed in NMR interferometry in 
a three-level system.^^ A superconducting phase qudit 
has been employed as an effective four-level system to 
show the symmetry of spinors under 27r-rotations, which 
can be interpreted as a geometric phase. ^ It has also 
been proposed to use three-level systems to detect non- 
Abelian geometric phases^^ and schemes to perform non- 
adiabatic holonomic quantum computation have been 
studied theoretically.^ 

Here, we present measurements of the effect of higher 
energy levels on the controlled accumulation of geomet- 
ric phases in a superconducting qubit system. In contrast 
to previous measurements of the geometric phase in su- 
perconducting circuits, where a Cooper-pair-box (CPB) 
qubit was used,^ we employ a transmon-type qubit. It 
is characterised by an increased ratio Ej/Ec of Joseph- 
son energy to charging energy of the qubit, a smaller 
charge dispersion^^ and a reduced anharmonicity, typi- 
cally only a few hundred MHz. As a consequence, the 
|1) ^ 1 2) transition of frequency CJ12 is in close vicinity 
to the |0) ^ |1) transition of frequency cJoi [Fig. [ija)]. 
Therefore, whereas the CPB qubit can safely be ap- 
proximated as a two- level system, the transmon cannot. 
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FIG. 1. (Color online) (a) Schematic energy level diagram 
of the transmon, with the computational subspace spanned 
by |0) and |1), and the second excited state |2) with anhar- 
monicity 0^2 < 0. (b) Lumped element circuit diagram of the 
sample and the measurement setup (see text for details), (c) 
Optical microscope image of the sample with two transmons 
coupled to a coplanar waveguide resonator with individual ca- 
pacitively coupled microwave gate lines, (d) Close-up of the 
transmon used in the experiments. 



Higher transmon levels affect the qubit dynamics, '^^'^ 
but can also be employed as a resource for quantum 
gates to improve single-shot readout, or to imple- 
ment a single-photon router. Quantum optical experi- 
ments involving three levels have been carried outP^H^ 
and their controlled preparation and tomography has 
been demonstrated."^ 

Here, we find that the difference in the level struc- 
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ture between the transmon and an effective two-level sys- 
tem also significantly affects the geometric phase. Time- 
independent perturbation theory is used to successfully 
model the geometric phase of a transmon. Then, the 
dependence of the geometric phase on evolution time in 
both the adiabatic and non-adiabatic regime is analysed. 
Non-adiabatic corrections to the geometric phase are ob- 
served and explained by simulating unitary dynamics. 

For manipulation and readout, the transmon is cou- 
pled to a coplanar transmission line resonator of quality 
factor Q = 2155 via a gate capacitance Cg. The lumped- 
element circuit diagram of the measurement setup, as 
well as an optical microscope image of the sample, is 
shown in Fig. [TJb-d). The sample is mounted in a 
dilution refrigerator and operated at a temperature of 
20 mK. The transition frequency of the transmon is tun- 
able by externally applied magnetic flux <l>, generated us- 
ing superconducting coils mounted underneath the sam- 
ple. The transmon state is manipulated using resonant 
and off-resonant microwave tones (of frequency cjqi and 
cjoi — A, respectively) created by AC modulation of an in- 
phase/quadrature- mixer. This provides individual con- 
trol of both quadratures {x and y) of the driving mi- 
crowave signal, which couples capacitively to the trans- 
mon via a local microwave gate line. 

From spectroscopic measurements, we have deter- 
mined a maximum Josephson energy Ej/h = 13.96 GHz 
and a charging energy Ec/h = 0.36 GHz, which corre- 
sponds to a maximum transition frequency cJoi,max/27r = 
5.95 GHz, and a coupling strength ^/27r = 360 MHz to 
the fundamental mode of the resonator. To reduce the 
Purcell effect and optimise coherence properties, cc;oi,max 
was designed to lie below the fundamental mode of the 
resonator.^ The experiment is carried out in the dis- 
persive regime, where the transmon, biased at co'oi,max, 
is non-resonantly coupled to the resonator (at frequency 
uOr/^Ti = 6.662 GHz with the qubit in the ground state) 
and can be read out via a state-dependent change in the 
microwave tone of frequency Ur transmitted through the 
resonator. Since Ej/Ec = 39 1 and the anhar- 
monicity is 0^2 /27r = —423 MHz, our sample is operated 
well within the transmon regime. Charge dispersion is 
expected to amount to about 100 kHz for cjqi and about 
2.9 MHz for CJ12. We have measured an energy relaxation 
time Ti = 0.84 /is, a phase coherence time T2* = 1.03 /is 
and a spin-echo phase coherence time Tf^^^ = 1.11 /is. 

We consider the Hamiltonian of a driven n-level qubit 
in a frame corotating at the drive frequency uJd^ 



i7((^) = n^(jA + a,)|j)(j| 
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where aj is the anharmonicity defined through the energy 
of the jth qubit energy level ujoj = joJoi + (with = 
Q^i = 0), is the strength of the drive expressed in units 
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FIG. 2. (Color online) (a) Slcetch of the microwave pulse se- 
quence used in a geometric phase measurement, consisting of 
a series of resonant (cjoi) and off-resonant (cjoi — A) pulses 
applied either on the x-quadrature (blue solid line) or the y- 
quadrature (green dotted line). The geometric phase is gen- 
erated by the adiabatic evolution of the qubit between the 
resonant pulses, (b) Simulated adiabatic evolution (blue line) 
of the ground state subjected to the off-resonant drive along 
the path C~, visualised in the Hilbert space of a three-level 
system and on the Bloch sphere, the approximate two-level 
equivalent (see text for details), (c) Extracted phase 7 as a 
function of solid angle A. Shown is the experimental data for 
C ^ (circles), C (diamonds) and (triangles), as well 

as the geometric phase obtained with second-order pertur- 
bation theory (solid lines) and the prediction for a two-level 
system (dashed lines). The off-resonant pulses were applied 
with detuning A/27r = —35 MHz. (d) Extracted phase 7 as a 
function of detuning A. The experimental data for solid an- 
gles A ^ 7r/4, 37r/4, and 57r/4 [indicated by circles, triangles 
and diamonds, respectively, and also indicated by arrows in 
(c)] is shown alongside the geometric phase calculated using 
second-order perturbation theory (solid lines) and the predic- 
tion for a two-level system (dashed lines). 



of angular frequency and cp its phase. A = cjqi — ojd 
denotes the detuning between the frequency of the |0) ^ 
|1) transition and the drive frequency. The increase of the 
coupling strength cx y^j + 1 follows from the off-diagonal 
matrix elements of the charge operator. In the two- 
level approximation, we restrict the Hamiltonian to the 
computational subspace of the qubit spanned by |0) and 
|1>: 

H ^^{Q,^a^ + Qyay + Aa,) = ^a-B, (2) 
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with Qx = ^cos(f and Qy = l^siiK/?. This is pre- 
cisely the Hamiltonian of a spin-half particle in an ef- 
fective magnetic field B = {Qx^ ^)-^ By applying mi- 
crowave frequency drives, we are able to experimentally 
control all three parameters of this effective field, which 
is our parameter space, to guide the qubit along a cir- 
cular path C with constant detuning A. At first, the 
drive field is ramped up, tilting B so that it forms an 
angle = arctan(l]/A) with respect to the z-axis. Then, 
the phase of the drive is swept by 27r, causing B to ro- 
tate once around the 2;-axis, either clockwise or anticlock- 
wise. Finally, the drive is ramped back to zero. The time 
evolution of both quadratures of the drive is plotted in 
Fig.[2|a). During this sequence, the solid angle subtended 
by the path as seen from the origin is A = 27r(l — cosi9). 
We then repeat the measurement for different driving 
strengths ft, thereby changing the solid angle A. 

To determine the geometric phase experimentally, we 
employ an interferometric measurement [Fig. [2|a)]. The 
leading and trailing resonant 7r/2-pulses implement a 
Ramsey measurement, while the resonant spin-echo tt- 
pulse in the centre serves to cancel the dynamic phase. ^ ^ 
After preparing an equal superposition of the |0) and |1) 
states, the qubit traverses the path C~ and the relative 
phase 2(^(1 — jg) acquired between |0) and |1) comprises 
both a dynamic (7^^) and a geometric (7^) contribution. 
The spin-echo 7r-pulse then effectively flips the sign of the 
phase. As it traverses the second loop in the opposite di- 
rection, the qubit acquires the phase 2(7^^ + 7^) since 
the dynamic phase, unlike the geometric phase, is inde- 
pendent of the direction of evolution. Thus, after follow- 
ing the contours C ^, dynamic phase contributions can- 
cel out and the qubit state has acquired a phase 7 = Ajg 
which is purely geometric. Tracing out the contours in 
opposite direction, simply inverts the sign of the 

phase, while following the contours twice in the same di- 
rection, C~^~^ or results in zero phase and serves as 
a control experiment [Fig. [2|c)]. 

During the off-resonant pulse sequences C^, the drive 
ft is strong (corresponding to induced Rabi-frequencies 
< 110 MHz) and therefore the higher levels of the qubit 
are populated as well. In order to visualise this popula- 
tion leakage, we consider the Hilbert space of a three- 
level system. Neglecting a global phase, any three- 
level state can be written as = e^^^ sin Pi cos/32|0) + 
giX2 sin/3i sin/32|l)+cos/3i|2), where xi^2 ^ [O7 ^tt] are the 
phases of the ground and first excited state, respectively, 
relative to the phase of the second excited state, and 
Pi,2 ^ [07 7r/2] parametrise the populations. Therefore, 
every state can be represented as a point in the product 
manifold of a torus and an octant of a unit sphere.!^ Ob- 
serving that /3i ^ 7r/2 while the ground state is subjected 
to the off-resonant drive C~ [Fig.[2|b)], we conclude that 
the instantaneous ground state leaves the computational 
subspace. The population of second excited state reaches 
up to ^ 12%, showing the necessity to consider the higher 
levels in our experiment. It is important to note that all 
resonant pulses effectively act in the subspace spanned 



by |0) and |1): in order to avoid exciting higher energy 
levels with the resonant pulses, we use an optimal con- 
trol technique known as DRAG.^ Furthermore, to en- 
sure that the second excited state |2) is depleted before 
the resonant pulses are applied, the off-resonant drive is 
adiabatically ramped down. 

After qubit manipulation, which takes approximately 
700 ns, the population = {l — crz)/^ of the first excited 
state is extracted by a dispersive readout.!^ The phase 7 
the qubit has acquired during evolution is reconstructed 
with state tomography (ST). The second 7r/2-pulse of the 
Ramsey sequence rotates the qubit about either the x or 
y axis, and serves as tomography pulse. In the absence of 
decoherence, the phase 7 is given by arctan ((cr^)/(craj)) 
with {(Jx) = cos 7, {(jy) = sin 7 and (a^) = 0. The same 
expression approximates 7 well even in the presence of 
decoherence, which reduces the size of the Bloch vector 
(da^, (Ty^az) to 0.47 in our experiments, while keeping the 
ratio (cTy) / {(Jx) constant. Therefore, the geometric phase 
remains unaltered by the decoherence in our experimen- 
tal setting. 

In keeping with Berry's predictions for a two-level sys- 
tem, 7 = 2A, we measure a phase 7 which is approx- 
imately twice the solid angle A subtended by the path 
[Fig. [2|c))]. However, the data clearly shows that 7 in- 
creasingly deviates from Berry's predictions as A (and 
therefore also the drive Q) increases: the measured ge- 
ometric phase is observed to be up to 15% larger than 
expected. 

These deviations can be explained by the presence 
of higher transmon levels. Defining the operator N = 
Yl^=oj\j){j\^ the Hamiltonian in Eq. (jlj) can be rewrit- 
ten as H{(p) = e-^^^i^(0)e^^^. It follows that, given an 
eigenvector |<l>(0)) of i^(0), \^{ip)) = e-^^^|^(0)) is an 
eigenvector of H{(p). For the circular path C described 
above, the geometric phase 7<^(o)_^cquired by the eigen- 
vector \^{0)) is then found to 



7<^>(o) 



-'f 

Jo 



(3) 



.(0) 



From Eq. (31), one indeed recovers the expression 7^^ 
7r(l ± cos'iTj for the geometric phase of a two- level sys- 
tem, where the sign ± corresponds to the positive and 
negative eigenvalue of i^(0), respectively. To compute 
the geometric phase for a multi-level system, we divide 
the Hamiltonian H{0) = Hq + V into a block-diagonal 
part Ho coupling the lowest two transmon levels, and a 
perturbative part V coupling the remaining levels: 



i7o = n^(jA + a,)|j)01 



(|l)(0|+h.c.). 



J=0 



^=^E(V7Tib' + i)01 + h.c.). 

Substituting the expansion of the eigenvectors in V, 
= + 1^5'^) + + into Eq. (g, and 
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FIG. 3. (Color online) (a) Measured phase 7 as function of 
phase sweep time r for the solid angle A = 7v/4 at detuning 
A/27r = —45 MHz. The dashed line is the phase obtained by 
numerically calculating the unitary time evolution of a four 
level transmon using the Schrodinger equation. The adia- 
baticity parameters a for a given r are indicated on the up- 
per horizontal axis, (b) Fidelity of the geometric phase gates 
shown in (a) as a function of r. 



retaining terms up to second order, the second excited 
state is found to contribute 



A7± 



: 7± — 7^^ = nk sin^ i9 

2A:(1±cos^) + (2A:t(3A:- 



(4) 



-2)cosi?) sin^i? 



(/cqp (3^ + 2) cos 



to the geometric phase 7±, where k = A/a2. In the ex- 
periment we measure the quantity A7 = 2(A7_ — A7+), 
see Fig. [SJc). We also note that A7 = 27rk sin^ i^/ cos 1? + 
0{k'^) vanishes for large a2 as expected. 

Since in the experiment /c ~ 1/8 is small and 1^1/^2 1 <C 



1, the expansion coefficients of 1^^^^) and l^Y^) are small, 
and perturbative treatment is justified. To verify the va- 
lidity of the perturbative treatment, we simulated the 
qubit evolution for the pulse sequence in Fig. |2|a), re- 
taining four energy levels in a quantum master equation 
simulation, thereby taking into qubit population decay, 
decoherence and non-adiabatic effects arising from finite 
evolution time.'^ Also, the Hamiltonian in Eq. ([l]) with 
n = 4 was numerically diagonalised to compute the ge- 
ometric phase in the limit of perfect adiabaticity. We 
found that the perturbatively computed geometric phase 
7 using Eq. Q differs from both simulations and nu- 
merical results by less than 2%. Furthermore, we have 
measured 7 for a range of detunings A and have found 
good agreement with the geometric phase computed us- 
ing perturbation theory [Fig. [2|d)]. 



We have verified that the geometric phase does not 
depend on the dispersive coupling of the transmon to 
the resonator by tuning cjqi such that S = ujr — ojqi = 
1.58 GHz and comparing the data to the case shown in 
Fig. ^c), where S = 0.71 GHz. 

Finally, the transition from the adiabatic regime to the 
non-adiabatic regime was examined by measuring the ac- 
quired phase 7 at fixed solid angle A = 7r/4, changing 
only the phase sweep time r [Fig. [2|a)] but keeping the 
duration of the total sequence constant. The data in 
Fig. [3]^ a) shows that 7 is constant, i.e., independent of 
evolution time for r larger than about 50 ns, or equiv- 
alently, for adiabaticity parameters a = (psin('i9)/|B| < 
0.1. The measured 7 also agrees well with the result ob- 
tained by numerically solving the Schrodinger equation. 
For values of a > 0.1, when entering the non-adiabatic 
regime, we observe that 7 oscillates and varies by more 
than 50%. In this regime, 7 is a combination of dynamic 
and geometric phase because the spin-echo technique fails 
to cancel the dynamic phase: for non-adiabatic evolution, 
the state after the spin-echo 7r-pulse does not necessar- 
ily correspond to the initial state with |0) and |1) inter- 
changed. 

In the context of quantum information processing, the 
manipulation sequence could serve as a single qubit phase 
gate. Its performance can be assessed by computing the 
fidelity F = tryp^^T?^, where p is the experimen- 
tal density matrix processed with maximum likelihoocPl 
and (J is the expected density matrix for perfectly adi- 
abatic evolution. We find that the fidelity of the gate 
averages F = 90% in the adiabatic regime [Fig. [sj^b)]. 
There, about 8% of the loss in fidelity can be attributed 
to qubit decay, whereas inaccuracies in qubit preparation 
and qubit dephasing account for the remaining 2%. In 
the non-adiabatic regime a significant decrease in fidelity 
is observed, as expected. 

In conclusion, we have measured the geometric phase 
in a multi-level system with small anharmonicity and ob- 
served that the two-level approximation breaks down for 
strong drives, as evidenced in our experiment by devi- 
ations of the geometric phase from the expected linear 
dependance on solid angle. We have modelled the contri- 
butions from the second excited state to the adiabatic ge- 
ometric phase using time-independent perturbation the- 
ory. By examining Berry's phase in the adiabatic limit 
and going to the non-adiabatic regime, we have shown 
that it is independent of evolution time for adiabaticity 
parameters ^0.1. The phase in the non-adiabatic limit 
could potentially inherit some of the adiabatic phase's 
noise resilience, suggesting further experimental tests on 
how it is affected by noise in the control parameters. 
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